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ABSTRACT 

Gas-dynamical studies using the Hubble Space Telescope are an integral component for future progress in the 
search for massive black holes in galactic nuclei. Here we present an extensive set of gas rotation curves obtained 
with the Space Telescope Imaging Spectrograph for the central regions of 23 disk galaxies. We find that the 
bulges of randomly selected, nearby spiral and SO galaxies generally do not contain well-defined gaseous disks. 
Only 15%-20% of disk galaxies have regular, symmetric velocity fields useful for dynamical analysis. Through 
comparison of the kinematics with Hubble Space Telescope images of the nuclear regions, we show that the 
probability of success can be significantly boosted by preselecting objects whose central dust lanes follow a well- 
ordered, circularly symmetric pattern. The dust morphology can be ascertained efficiently by visual inspection of 
unsharp-masked images. 

Subject headings: galaxies: bulges — galaxies: ISM — galaxies: kinematics and dynamics — galaxies: nuclei 
— galaxies: spiral — galaxies: structure 



\. INTRODUCTION 

The installation of the Space Telescope Imaging Spectro- 
graph (STIS; Leitherer et al. 2001) on the Hubble Space Tele- 
scope (HST) ushered in a new era in the search for massive 
black holes (BHs) in galactic nuclei. By now, sufficiently large 
numbers of BH candidates are known that the emphasis has 
shifted from detections of individual objects to demographic 
studies using statistical samples (e.g., Richstone et al. 1998; 
Ho 1999; Gebhardt et al. 2000b; FeiTarese & Merritt 2000). 
Future progress will focus on two fronts, first on sharpening the 
analysis techniques to improve the accuracy of the mass mea- 
surements, and second on substantially enlarging the existing 
samples to cover galaxies across the Hubble sequence. 

BH searches with HST use either stars or gas to trace the cen- 
tral potential. These two approaches are complementary, each 
with its merits and limitations. Modern, fully general, stellar- 
dynamical models are powerful because they give information 
not only on masses but also on the orbital structure of the galaxy 
(e.g., van der Mai-el et al. 1998; Cretton & van den Bosch 1999; 
Gebhardt et al. 2000a). One drawback is that both the data and 
the computational requirements are expensive. Another is that 
the modeling can be very challenging for complicated dynam- 
ical structures; the double nucleus of M31 serves as a good il- 
lustration (e.g., Kormendy & Bender 1999; Bacon et al. 2001). 
This technique also cannot easily cope with dusty galaxies, as 
a consequence of which most stellar-dynamical work with HST 
has been biased toward early-type (mostly elliptical) galaxies. 

Gas kinematics in principle offer a more straightforward 
measure of the central mass, provided that the gas participates 



in Keplerian rotation in a disklike configuration. The observa- 
tions are less demanding because optical emission lines from 
ionized gas generally have much higher equivalent widths than 
stellar absorption lines. Moreover, as demonstrated later in this 
paper, optical nebular emission invariably traces nuclear dust 
features. This technique, therefore, can be applied to galaxies 
inaccessible to stellar observations. Unlike stars, however, gas 
behaves as a collisional fluid and responds to nongravitational 
perturbations such as shocks, radiation pressure, or magnetic 
fields. Despite these potential complications, gas-dynamical 
models have been applied to a number of elliptical galaxies 
with circumnuclear gas disks (Harms et al. 1994; Ferrarese, 
Ford, & Jaffe 1996; Macchetto et al. 1997; van der Marel & 
van den Bosch 1998; Bower et al. 1998; Ferrarese & Ford 1999; 
Verdoes Kleijn et al. 2000). 

Our group has concentrated on using STIS to apply the gas- 
dynamical technique to the bulges of spiral and SO galaxies 
(Shields et al. 2000; Barth et al. 2001; Sai'zi et al. 2001, 2002). 
Our work has revealed a subtlety that has not been widely ap- 
preciated in previous studies: the majority (^^80%) of randomly 
chosen disk galaxies do not have well-behaved velocity fields 
suitable for gas-dynamical modeling at HST resolution. Based 
on this experience, we have developed an efficient strategy to 
preselect galaxies that significantly increases the chance of suc- 
cess for future STIS programs. This paper describes our tech- 
nique, which is based on the use of the dust morphology to 
predict the velocity field. 
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Fig. la. — Dust morphology and gas kinematics for NGC 278, 2787, 3351, 3368, 3489, and 3982. For each galaxy, the top panel shows the unsharp-masked, 
central 9f' 1 x 9^' 1 region of either a V-band WFPC2 image or the F28X50LP STIS acquisition image. (For NGC 3368, the image is a 20" x 20" section of the WF2 
chip.) The horizontal bar on the upper right comer gives the scale of 1". The rectangular region superposed on the nucleus outlines the position of the slit. The 
bottom panel shows the Ha velocity curve of the central ~0"5-l"0, where the rest-frame velocity of the galaxy was determined by centering the rotation curve on 
the peak of the stellar continuum. The velocity curve is omitted for galaxies lacking extended emission. 
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Fig. lb— Dust morphology and gas kinematics for NGC 3992, 4143, 4203, 4245, 4314, and 4321. As in Figure la. 
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Fig. Ic— Dust morphology and gas kinematics for NGC 4380, 4435, 4450, 4459, 4477, and 4501. As in Figure la. 
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Fig. Id — Dust morphology and gas kinematics for NGC 4548, 4596, 4698, 4800, and 5055. As in Figure la. 
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TABLE 1 
The Sample 



Galaxy Name 


Hubble Type 


T 


P.A. 


i 


D 








(deg.) 


(deg. 


) (Mpc) 


(1) 


(2) 


(3) 


(4) 


(5) ' 


(6) 


NGC 278 


SABb 


3.0 






11.8 


NGC 2787 


SBO+ 


-1.0 


117 


51 


13.0 


NGC 3351 


SBb 


3.0 


13 


48 


8.1 


NGC 3368 


SABab 


2.0 


5 


47 


8.1 


NGC 3489 


SABO+ 


-1.0 


70 


56 


6.4 


NGC 3982 


SABb: 


3.0 




30 


17.0 


NGC 3992 


SBbc 


4.0 


68 


53 


17.0 


NGC 4143 


SABO 


-2.0 


144 


52 


17.0 


NGC 4203 


SABO-: 


-3.0 


10 


21 


9.7 


NGC 4245 


SBO/a: 


0.0 




41 


9.7 


NGC 4314 


SBa 


1.0 




27 


9.7 


NGC 4321 


SABbc 


4.0 


30 


32 


16.8 


NGC 4380 


Sb? 


3.0 


153 


58 


16.8 


NGC 4435 


SBO 


-2.0 


13 


44 


16.8 


NGC 4450 


Sab 


2.0 


175 


43 


16.8 


NGC 4459 


S0+ 


-1.0 


110 


41 


16.8 


A AT'? 

iN^jT^^ 44/ ; 


SBO? 


-2.0 


15 


0/1 
Z4 


1 Q 


NGC 4501 


Sb 


3.0 


140 


59 


16.8 


NGC 4548 


SBb 


3.0 


150 


38 


16.8 


NGC 4596 


SBO+ 


-1.0 


135 


43 


16.8 


NGC 4698 


Sab 


2.0 


170 


53 


16.8 


NGC 4800 


Sb 


3.0 


25 


43 


15.2 


NGC 5055 


Sbc 


4.0 


105 


56 


7.2 


Note.— Col. (1] 


1 Galaxy name. Col. (2) 


Hubble type 


. Col 


. (3) Morpho- 


logical type index. 


Col. (4) Position 


angle of the major 


axis 


of the galaxy. 



Col. (5) Inclination angle of the galaxy. Col. (6) Adopted distance, based 
on Ho = 75 km s"^ Mpc"^ Data taken from Ho et al. 1997. 



2. THE DATA 

The data used in this paper derive from our Cycle 7 STIS pro- 
gram GO-7361 . Since a full description of the observations and 
data reductions will be given elsewhere (Rix et al. 2002; see 
also Sarzi et al. 2001, 2002), here we only mention a few key 
points pertinent to the present discussion. Our sample was de- 
signed to address a diverse set of scientific goals, among them 
measurement of BH masses in bulges. It comprises nearly all 
northern, bright SO-Sb galaxies within 17 Mpc that are known 
to have nuclear emission in the spectroscopic catalog of Ho, 
Filippenko, & Sargent (1997). We avoided galaxies later than 
Sb because their bulges are expected to be faint or absent, and 
because their higher dust content may compromise our ability 
to center the narrow slit on the nucleus. To further ensure that 
the nucleus wiU be easily accessible, we eUminated galaxies 
more inclined than 60°. The visibility of the nucleus was also 
confirmed by inspection of WFPC2 images, which were avail- 
able for most of the galaxies. The final sample contained 24 
galaxies, of which 23 (Table 1) were successfully observed; the 
acquisition of NGC 4138 failed because a nearby, bright star 
was mistaken for the nucleus. 

Following a brief target-acquisition image using the 
F28X50LP long-pass filter, the primary kinematics data for 
each galaxy consist of a single-orbit integration with the 
G750M grating (AA6300-6870 A; FWHM AA w 2.2 A) using 
the 0"2x52" slit. To enable flexibility in scheduling, we did 
not choose specific position angles for the slit. 



3. DUST MORPHOLOGY AND VELOCITY FIELD 

We detected Ha and [N IIJ AA6548, 6583 emission in all flie 
objects, but the emission was symmetrically extended in the 
central 0."5-l."0 in only ~75% (17/23)ofthe sample. Forfliese 
17 objects, "rotation curves" were derived from the radial ve- 
locities of the gas as a function of position along the sUt using 
the Ha and [N II] emission lines. The velocities of these fea- 
tures were derived using SPECFIT (Kriss 1994) within IRAF^, 
assuming a common velocity, a Gaussian profile, and a fixed 
[N II] doublet ratio of 3 : 1 . The results are shown in Figure 1 . A 
striking result is the rarity of galaxies with well-behaved veloc- 
ity fields. Among the 17 galaxies with spatially resolved veloc- 
ity profiles, only four (NGC 2787, 4203, 4459, 4596), or 25%, 
show ordered rotation curves amenable to dynamical analysis 
(Sarzi et al. 2001). The rest of the velocity curves look disor- 
ganized and difficult to interpret (Fig. 1). In terms of using the 
data for BH searches, therefore, the success rate is disappoint- 
ingly low. 

The greyscale panels in Figure 1 display the unsharp-masked 
optical continuum image of each galaxy. Most of the data are 
archival broad-band (roughly V) images acquired with the PC 
detector on WFPC2; two are our own STIS F28X50LP target- 
acquisition images. To accentuate high-spatial frequency fea- 
tures such as star clusters and dust lanes, we subtracted from 
each image a version of itself smoothed by a Gaussian with a 
= 4 pixels (^0"2). The nucleus, by design, is clearly visible in 
each image. As was known from previous //5r imaging studies 
of disk galaxies (Phillips et al. 1996; CaroUo et al. 1997), dust 
is ubiquitous. Every galaxy in our sample shows dust lanes at 
some level within a central radius of ~5", or ~400 pc for a 
typical distance of 17 Mpc. However, unlike the case of early- 
type (E and SO) galaxies (van Dokkum & Franx 1995; Tomita 
et al. 2000; Tran et al. 2001), the dust distributions do not 
resemble simple disks with sharply delineated, outer edges. In- 
stead, the dust morphology can be best described as nuclear spi- 
rals (e.g., NGC 3351, NGC 3982), concentric rings (e.g., NGC 
2787, NGC 4596), or simply chaotic (e.g., NGC 3368, NGC 
4450). The one object with a fairly well-defined nuclear disk, 
NGC 4435, turns out to be an SO galaxy. 

Closer examination of the data reveals that the degree of "or- 
derliness" of the gas kinematics evidently correlates with the 
morphology of the dust lanes. All four galaxies with well- 
behaved velocity fields (NGC 2787, 4203, 4459, 4596) have 
relatively smooth, circularly symmetric dust rings. The con- 
verse seems to hold for most, but not all, objects. For example, 
die dust patterns in NGC 4245 and NGC 4800 resemble those 
in the above-mentioned four objects, but neither has a regular, 
symmetric velocity field in ionized gas. The case of NGC 4435 
deserves special remark. Although this galaxy possesses a clean 
nuclear disk, we were unable to obtain a usable symmetric rota- 
tion curve because the slit fell near the minor axis of its highly 
inclined disk (Fig. Ic). The velocity curve of the far side of the 
disk indeed looks regular, but that of the near side is obstructed 
by a prominent dust lane. 

Thus, in terms of planning for HST BH searches using gas 
kinematics, it appears that the morphology of the dust lanes can 
be used as an efficient and relatively reliable predictor of the 
gas kinematics. Albeit somewhat subjective, and by no means 
perfectly foolproof, the strategy we advocate is nonetheless 
a highly effective tool for cuUing promising targets for STIS 



' IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under 
cooperative agreement with the National Science Foundation. 
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spectroscopy. It is especially attractive in view of the extensive 
database of optical images of nearby galaxies available in the 
HST archives and the straightforward manner in which unsharp 
masking can be performed on them. It would be fruitful to re- 
visit these issues in the future when we have available images 
and spectra with resolutions higher than that currently possible 
with HST. 

4. SUMMARY 

Searches for massive BHs in bulges of disk galaxies gener- 
ally must rely on gas-dynamical methods because these objects 
are dusty. Our experience with STIS data indicates that the suc- 
cess rate will be low in randomly chosen samples of spiral and 
SO galaxies because well-defined nuclear gas disks with regu- 
lar, symmetric velocity fields are rare in these objects. We show 
that the degree of regularity in the velocity field correlates with 
the morphology of the dust lanes, which can be efficiently as- 



certained through unsharp masking of optical continuum HST 
images. Objects with smooth, circularly symmetric dust lanes 
tend to have velocity fields that are amenable to dynamical anal- 
ysis. By contrast, galaxies with obviously disorganized dust 
morphology will almost never yield easily interpretable veloc- 
ity fields. We draw attention to this empirical correlation as 
a powerful tool to select targets for future HST programs to 
search for massive BHs using gas kinematics. 



This work was supported financially through NASA grant 
NAG 5-3556, and by grants in support of HST programs GO- 
7361 and GO-7403, awarded by STScI, which is operated by 
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